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Abstract. We study the dynamics of polymer chains in a bath of self-propelled 
particles (SPP) by extensive Langevin dynamics simulations in a two dimensional 
system. Specifically, we analyse the polymer looping properties versus the SPP activity 
and investigate how the presence of the active particles alters the chain conformational 
statistics. We find that SPPs tend to extend flexible polymer chains while they rather 
compactify stiffer semiflexible polymers, in agreement with previous results. Here we 
show that larger activities of SPPs yield a higher effective temperature of the bath 
and thus facilitate looping kinetics of a passive polymer chain. We explicitly compute 
the looping probability and looping time in a wide range of the model parameters. 
We also analyse the motion of a monomeric tracer particle and the polymer’s centre 
of mass in the presence of the active particles in terms of the time averaged mean 
squared displacement, revealing a giant diffusivity enhancement for the polymer chain 
via SPP pooling. Our results are applicable to rationalising the dimensions and looping 
kinetics of biopolymers at constantly fluctuating and often actively driven conditions 
inside biological cells or suspensions of active colloidal particles or bacteria cells. 


1. Introduction 

Active motion is a necessary prerequisite for living systems to maintain vital processes, 
including materials transport inside cells and the foraging dynamics of mobile organisms 
0 0. The length scales associated with active motion processes span several orders 
of magnitude and range from the nanoscopic motion of cellular molecular motors 
[3] essential to move larger cargo in the crowded environment of cells |1], over the 
microscopic motion of bacteria cells and micro-swimmers 016], to the macroscopic 
motion patterns of higher animals and humans [7]. In particular, artihcial Janus colloids 
are propelled by diffusiophoretic or thermophoretic forces mu [To] . Active motion 
enhances the speed and precision of signalling and cargo transport in biological cells 
m and allows efficient search of sparse targets for large organisms [12]. 
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A somewhat different question is how passive particles are influenced by an active 
environment. Tracking the motion of tracer particles immersed in baths of active 
bacteria and swimming eukaryotic cells na one typically observes an enhanced 

effective tracer mobility, and the active environment may lead to exponential tails of 
the displacement distribution [151 IB]- Passive particles may also become enslaved to 
the motion of motor-cargo complexes due to cytoplasmic drag na. When micron sized 
colloids are immersed in baths with smaller particles, short ranged attractive depletion 
forces of entropic origin are observed [IH]. However, the same colloids in a bath of 
self-propelled particles (SPPs) may experience long ranged attractive or repulsive forces 
depending on the SPP characteristics ng. By tuning the concentration of SPPs the 
forces between two plates can be controlled [20] , 

Here we want to focus on the properties of polymer chains in an active liquid. 
It is known that when a polymer chain is immersed into an SPP bath its extension 
changes non-monotonically with the activity Fa of active particles due to competing 
effects of active forces and chain elasticity [211122]. We study here the extent to which 
the activity of SPPs alters the internal motion of a polymer chain, specifically, how 
its end loop formation kinetics becomes affected. Polymer looping or reactions of the 
chain ends is a fundamental process governing numerous biological functions [23l | 23 |. 
Protein mediated DNA looping, for instance, is involved in gene regulatory processes 
[2511261127], or DNA and RNA constructs may be used as molecular beacon sensors |28j . 

We quantify the behaviour of a polymer chain immersed in a bath of SPPs in two 
dimensions, see the snapshot of our systems in hgure [H We hnd that the activity Fa 
of SPPs differently affects the chain conformations depending on the chain bending 
stiffness k. Concurrently, SPPs significantly enhance the looping kinetics as well as 
give rise to a giant diffusivity of the centre of mass motion of the chain due to SPP 
pooling in typical parachute like chain conhgurations. We analyse the diffusion of a 
monomeric tracer particle, which shows superdiffusive motion on short time scales and 
Brownian behaviour with enhanced diffusivity in the long time limit. In references 
[ZD [22] a similar system was considered, the main focus being on equilibrium polymer 
properties such as the gyration radius of the chain. Below we systematically analyse 
dynamic properties of the polymer chain. Our results demonstrate that the equilibrium 
and dynamic properties of polymers in an SPP bath are to be considered on the same 
footing. 

This paper is organised as follows. We introduce our model and the simulations 
methods in section [D In section 13.11 we examine the equilibrium properties of the 
polymer chain. Section 13.21 presents the main results regarding the polymer looping 
properties. In section l373l we study the dynamical effects of SPPs on the tracer diffusion, 
in order to understand its implications on the enhancement of the polymer looping 
kinetics. We summarise our results and discuss their possible applications in section jH 
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To study the dynamics of a polymer chain immersed in a bath of SPPs, we employ 
coarse grained computer simulations. The polymer chain is modelled as a bead spring 
chain consisting of n monomers of diameter a, connected by harmonic springs with the 
potential 




i=2 


ri-i| 


h 


( 1 ) 


Here k is the spring constant and /q is the equilibrium bond length. The self avoidance 
of the chain monomers is modelled by the repulsive part of the Lennard-Jones (LJ) 
potential (the so called Weeks-Chandler-Andersen or WCA potential), 

12 /cr\6l 


Uu{r) = 4e 


a 


+ C'(l"cut)- 


( 2 ) 


for r < Tcut, where rcut is a cutoff length. Moreover C(rcut) is a constant that ensures 
that Ui^j{r) = 0 for separations r > r^ut, and r is the inter-monomer distance. The 
potential strength is denoted by e. With the standard choice rcut = for the cutoff 
length the potential is purely repulsive. In what follows, we measure the length in units 
of a and the energy in units of the thermal energy ksT, where ks is the Boltzmann 
constant and T is the absolute temperature. Below we set the model parameters to 
a = 1, Iq = 1.12, k = 10^, and e = 1. 

The bending energy of the chain is given by 


11 ± 2 
^ i=2 

where k is the bending stiffness. For a given value of k, the chain persistence length 
is Ip ~ 2k/q/(/cbT) in two dimensions. The end monomers are subject to short ranged 
attractive interactions with energy e^, mimicking the biologically relevant situation that 
closed structures are energetically prohtable, as known for specihc DNA looping [23] or 
closed single stranded DNA (hairpins) [29]. We include the attractive interactions via 
the LJ potential in Eq. ([2]) but with a larger cutoff distance and attraction strength e^, 
namely 

U^ir) = t/Lj(r, e^) 


and rcut = 2cr. The effects of the end-to-end stickiness on the looping properties were 
considered by us recently [30]. In what follows we set Cs = bk-sT. 

The dynamics of the position rj(t) of the Ah chain monomer is described by the 
Langevin equation 

= -V|U + f/u(r) + f/J - 7^ + «.(*). (4) 

Here m is the monomer mass, 7 is its friction coefficient coupled to the diffusivity via 
D = /cbT/ 7 , (5) 
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Figure 1. Typical conformation of a polymer (blue chain) in a bath of active particles 
(red discs) in two dimensions. Parameters: the polymerisation degree of the chain is 
n=32, the packing fraction of SPPs is (/)=0.1, and the particle activity is Fq= 40 (see 
text for details). Video files for different chain lengths and SPP activities are available 
in the Supplementary Material. 


and ^i{t) represents a Gaussian white noise source of zero mean with autocorrelation 
{^i(t) ■ ^ii(t')) = ‘FykBT5i^i'5{t — t'), where 5*,*/ is the Kronecker delta symbol. 

The SPPs are modelled as disks of diameter a moving under the action of a constant 
force along a predehned orientation vector 


iij = {cos(6'j), sin(6ij)} . 


( 6 ) 


SPPs interact with each other as well as with polymer monomers via the WCA potential 
(j^. and the position of each SPP is governed by the Langevin equation [22] 

(PrAf) __ , , ^ , , dr, ^ , 

^ + Faiijit) - 7^ + ^j{t). (7) 

Here Fa is the active force amplitude, which is directly related to the SPP propulsion 
strength: it can be expressed in terms of the Peclet number Pe and the particle velocity 
V in terms of 


Pe = 


va 

15 


Fga 

ksT' 


( 8 ) 


The orientation 6j of the velocity of the jth SPP is changing as function of time according 
to the standard stochastic equation 


6 = \/2Dj. X ^r{t). 


(9) 
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Figure 2. Probability distribution of the polymer end-to-end distance for different 
SPP activities Fa and varying chain bending stiffness: (A) k = 0, (B) k = 3, and (C) 
K = 9. In (D) the looping probability Pi is plotted as function of Fa for n = 32 and 
= 0.05. 


Here fr is the Gaussian white noise associated with the rotational diffusivity Dy. which 
satishes the relation Dr = 3D/a^ (see, for instance, reference [22] )• Passive particles 
correspond to Fa = 0 , the situation studied in the context of polymer looping with 
macromolecular crowding in reference PI. 

In our simulations we use periodic boundary conditions for a square box of area 
where, depending on the length of the simulated chain, L varies from 60 to 80. The 
packing fraction of SPPs is dehned as </> = NcrAcr/L'^, where Ncr is the number of SPPs 
and Act = Ti{a/2)‘^ is the area of a single SPP. We consider rather dilute SPP systems 
with 0 < 0.1. For both chain monomers and active particles we choose the unit mass 
m = 1 and a relatively large friction of 7 = 5 to ensure a quick momentum relaxation. 
The time scale in the system is set by the elementary time to = cx^m/(/ cbT) [M]. We 
implement the Verlet velocity algorithm [3T] to simulate equation (|1]) and equation (17|) . 
The integration time step is At = 0 . 002 , and we typically run ~ lo® -^ steps to compute 
the quantities of interest. 

Generally the activity of SPPs may vary, or some particles in the bath may be 
completely inactive. To account for this fact, in a part of our study below we consider 
mixtures of active and inactive particles with respective fractions (fa and 0*. All these 
particles have the same mass and radius in the simulations. 
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Figure 3. Distribution of the gyration radius for varying Fa (in each panel) and 
bending stiffness (A) k = 0, (B) k = 3, and (C) k = 9. (D) Average gyration radius 
as function of the SPP activity. Parameters: n = 32 and f = 0.05. 


3. Main Results 
3.1. Polymer Dimensions 

We first consider the probability distribution function (PDF) of the end to end distance 
r of the chain as extracted from long time computer simulations, see hgure [2J4-C. In 
our simulations, due to the attraction between the end monomers the standard PDF 
of the polymer [27] acquires an additional peak around the minimum of the attraction 
potential at the end to end distance r ~ For a flexible chain with k = 0 (hgure 

|2jA) the chain gets more extended and the peak of the PDF shifts to larger distances 
when the activity Fa of SPPs increases. Conversely, for semihexible chains with a hnite 
value of the polymer stiffness k > 0, the peak is shifted to shorter distances (hgure 
[2B,C). These trends are similar to those of reference [19]. 

This is the main ehect of active particles on the static properties of passive polymer 
chains in solutions. Inspecting snapshots of the simulations (see also hgure H]) or the 
video hies in the Supplementary Material, one recognises that active particles ehect U 
or parachute like shapes of the polymer. Such parachute shapes are also observed for 
membranous red blood cells in cylindrical capillary hows and in blood vessels, see, e.g., 
reference [32]. For larger Fa values, when the SPP forces are much larger than the 
energetic scale for polymer bending, the PDF of the polymer end to end distance is 
nearly independent on the chain stihness k. 

In hgure [3] we also show the distribution of the radius of gyration of the 
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chain and its average valne (Rg)- For flexible chains the PDF of the gyration radius 
broadens towards larger values, causing the monotonic increase of {Rg)- In contrast, for 
semiflexible chains the gyration radius decreases for Fa < 20, above this value it only 
slightly increases, compare figure [3p. This behaviour is consistent with the results of 
references [211 [22]. 

3.2. Looping Probability and Looping Time 

The polymer end to end distance shows a highly erratic behaviour as function of time 
(see also the movies in the Supplementary Material for chains of different flexibility). 
The polymer ends tend to remain at short distances ~ due to their attractive energy, 
while longer end to end distances with r ~ rgq are favourable entropically at equilibrium. 
We compute the looping probability Pi and the looping time Ti from the time series of 
the polymer end to end distance r{t) generated in simulations. Similar to our recent 
studies [SO] [33] the looping probability Pi is defined as the fraction of time during which 
the end to end distance of the chain is shorter than r^. In that sense the critical distance 
Tc = 1.75cr separates the looped and un-looped states of the polymer, compare figure 
EK-C. 

Figure [2D shows the looping probability Pi as function of the activity Fa of 
the SPPs. For flexible chains the value of Pi decreases monotonically with Fa. 
Conversely, for semiflexible polymers the looping probability is non-monotonic in Fa. 
This observation indicates two competing effects of the active particles: on the one hand 
SPPs increase the effective chain flexibility resulting in higher Pi values. On the other 
hand, SPPs facilitate the unbinding of end monomers. We observe that, consistent with 
the shape of the end to end distance PDF at large activity Fa of SPPs in figures [2J4-C, 
for large Fa the looping probability is almost independent of k, as demonstrated by 

figure ED- 

The polymer looping time Ti is defined as the time interval within which the distance 
r reaches rgq for the first time and the time it gets shorter than the final distance 
rf = 1.2a, details are shown in figure 3 of reference [30]. While the distances and rj 
are mainly determined by the properties of the attractive potential of the end monomers, 
the value of rgq strongly varies with the chain length and the SPP activity Fa. From 
the PDFs of the end to end distances we first determine rgq for a given chain length na 
and particle activity and then use them to compute the looping time p. 

Although the effects of SPPs onto the spatial extension of the immersed polymers 
were considered previously [211 [22], their dynamic effects—particularly on the polymer 
end looping reaction—were not addressed in detail. In figure [Tj we present the polymer 
looping times as function of the particle activity Fa. In free space or for Fa = 0 
the looping takes much longer for stiffer chains because of the large bending energy 
required for a loop formation. As the SPP activity Fa increases the polymer looping 
time decreases, especially for stiff chains: we observe a reduction of Ti of more than two 
orders of magnitude, as evidenced in figure [H Figure E] shows that for large Fa values 
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Figure 4. Polymer looping time versus SPP activity Fa for polymer chains with n=32 
monomers. The simulation time of one point in this figure on a standard workstation 
is around 6 h. The error bars for Pi (figure [2p) and Ti were computed as the standard 
deviation of the mean via subdividing the time traces into ten subsets. 


the looping probability of flexible and semiflexible chains behaves quite similarly, and 
hgure m demonstrates similar trends for the looping times Ti of the polymer chains for 
large SPP activities. 

Up to now, we only considered chains of length n = 32 monomers. In hgure [5] we 
now show the looping times as a function of the chain length n. In free space (0 = 0), 
the looping time follows the scaling behaviour BOI 

( 11 ) 

with the Flory exponent v = 3/{d + 2) = 3/4 for a polymer in two dimensions. With 
= 0 (inactive crowders) the looping time increases for a given chain length mainly 
due to a decreasing monomer diffusivity in the medium [30]. In the presence of active 
crowders the polymer kinetics becomes facilitated, especially for longer chains, as shown 
by the red dots in hgure IH Interestingly, in the presence of active particles, the scaling 
exponent of Ti{n) decreases somewhat as compared to free space and passive crowders. 

In hgure [6] we show the polymer looping time versus the relative fraction of active 
particles, (pa/<P for the total crowding fraction p = pa + <Pi = 0.05. We observe that for 
small values pa/p the magnitude of Ti initially drops sharply, while the decrease of Ti for 
larger fractions of active particles is rather moderate. This indicates that the transition 
from the non-active to the active results in hgure |5] is rather non-uniform when active 
particles are added into the solution. 

3.3. Tracer Diffusion, Polymer Diffusion and Monomer Displacements 

To get a feeling for the ehects of SPPs on the dihusion of passive particles we now 
quantify the enhancement of the dihusive motion of a non-active tracer in a bath of 
SPPs. We track a particle with diameter a (same size as the monomers and SPPs) 
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Figure 5. Looping time versus chain length, plotted for the situation in free space 
as well as for conditions of crowded inactive (blue) and active (red symbols) particles. 
Parameters: (^=0.05 and k = 0. The asymptotic behaviour of equation (HU is shown 
as the solid line, the slopes 2.2 and 2.7 are also shown. 



Fraction of active particles 


( 12 ) 


Figure 6. Polymer looping time versus fraction fa/f of SPPs at total crowder fraction 
(/)=0.05. The chain length was n = 32 and the bending stiffness k = 3. 


for varying SPP activities Fa- From the time series of the tracer particle position 
r(t) = {x(t),y(t)} generated in our simulations we calculate the time averaged mean 
squared displacement (MSD) [SU [35] 


<5^(A) = 


T- A 


r -| 2 

x{t' + A) — x{t') dt', 


(13) 


where A is the so called lag time dehning the width of the averaging window shifted 
along the trajectory. Hereafter, the time averaged MSD is computed with respect to 
one dimension only. The additional mean 


1 ^ 
2=1 


( 14 ) 
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Time, t 

Figure 7. Time averaged MSD of a single tracer at varying SPP activities Fa 
plotted for </> = 0.05. Inset: effective tracer diffusivity, plotted in log-linear scale and 
normalised to the free space value. 


over an ensemble of N individual traces Sf{A) will be analysed below. 

We present the time averaged MSD of the tracer in £gure[7]for varying SPP activity 
Fa- The time averaged MSD grows faster than for Brownian motion (superdiffusion |34] ) 
only at very short times, A* < 2 ... 5, and then turns into the linear Brownian scaling, 
as expected. Since the momentum relaxation time, dehned as ~ m/y = 0.2, is shorter 
than the time scale A*, the extended super diffusion regime is likely due to the impact 
of active particles. 

We extract the diffusivity of the tracer particle through a linear £t to the long 
time behaviour of the time averaged MSD in the range A = lO^ '^. As shown in the 
inset of figure [7]in log-linear scale the diffusivity increases almost exponentially with Fa- 
This enhancement is the main reason of the dramatic facilitation of the polymer looping 
kinetics by highly active particles, as demonstrated in hgure 0] as function of the SPP 
activity Fa- This is one of the main conclusions of the current paper. 

Similarly, in figure [HJ4 we show the time averaged MSD of the polymer chain’s 
centre of mass for different SPP activities Fa- Comparison of the magnitudes of the 
time averaged MSDs shows that, as expected, the centre of mass diffusion of the entire 
chain is evidently much slower than that of a single tracer particle. In hgure [HB we 
compute the local scaling exponent of the MSD [SU 135] 


m) = 


d\og (5^(A) 


(15) 


(ilog(A) 

We observe that at short time scales the MSD increases superdiffusively with (3 > 1 and 
the anomalous diffusion exponent grows with increasing Fa values. At longer times the 
diffusion exponent decreases and around A ~ 10^ the motion becomes nearly Brownian, 
albeit with an enhanced diffusivity at higher Fa values. 

In the inset of hgure IHB we show the chain dihusivity of the centre of mass motion 
as function of Fa, normalised with respect to the value in free space. The enhancement 
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Figure 8. (A) Time averaged MSD of the polymer centre of mass motion, averaged 
over N = 5 trajectories, as a function of SPP activity. Inset: Effective diffusivity of the 
centre of mass motion, normalised to the free-space value. (B) Local scaling exponent 
/3(t) of the time averaged MSD. Parameters: n=32, k = 0, and (f) = 0.05. 


of the diffusivity of the polymer centre of mass is nearly two orders of magnitude, that 
is much higher than that of a monomeric tracer particle shown in the inset of hgure 
[71 This stronger enhancement is due to pooling of SPPs in the concave region of the 
parachute-shaped polymer chain, resulting in directed motion and faster diffusion of 
the polymer. This remarkable hnding of a giant diffusivity enhancement is our second 
major result. 

We also show the PDFs of the displacement of the polymer chain and of the tracer 
particle in panels A and B of hgure [9l respectively. Both for active and inactive crowders 
the PDFs exhibit Gaussian prohles, see the dotted hts in the hgure. In the presence 
of active particles, the width of the corresponding displacement PDF becomes wider, 
consistent with the enhanced dihusivity. This is particularly clear for the polymer centre 
of mass displacement shown in panel A of hgure [9l Interestingly, even at short times— 
when the time averaged MSD shows a superdihusive scaling—the PDFs remain nearly 
Gaussian. 
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Figure 9. PDFs of the centre of mass displacement of the polymer chain (A) and 
tracer particle (B), plotted for a set of diffusion times t and for two SPP activities at 
crowding fraction (/)=0.05. The dotted lines are Gaussian fits. 


4. Conclusions 

Actively driven systems are inherently out of equilibrium and exhibit peculiar 
behaviours, for instance, in the ratcheting of motors [36], the formation of living 
crystals [ST] |3H], phenomena of ordering [39], mesoscale turbulence phenomena 001 , or 
superfluidic behaviour may be observed in bacterial suspensions [H]. Even elementary 
laws of thermodynamics may no longer hold [12] SS] 01] . In that sense the behaviour of 
active liquid systems is as rich as that of active soft matter [15] . 

We studied the dynamics of a polymer chain in a bath of SPPs using Langevin 
dynamics simulations in two dimensions. We hrst considered the equilibrium behaviour 
of the gyration radius, the end to end distance distribution, and the looping probability 
of the chain as function of the activity of the SPPs. We found that a flexible polymer 
extends monotonically with the activity. In contrast, for a semiflexible chain—due 
to a competition of the the chain bending and active forces—the polymer size varies 
non-monotonically with the particle activity. For a larger activity of SPPs—when 
active forces become dominant over the chain elasticity effects—the extension of both 
semiflexible and flexible chains behaves quite similarly. 

SPPs also significantly impact the polymer kinetics, the focus of this study. Overall, 
due to the enhanced diffusivity of the chain monomers, the looping dynamics becomes 
faster. Especially for the case of stiffer chains the presence of SPPs renders the chains 
effectively softer, and the looping kinetics is dramatically enhanced in comparison to that 
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of flexible chains. Our results indicate that the activity of a cell medium, as mimicked 
above by active particles, may indeed facilitate DNA loop formation, effectively making 
the molecule more flexible. 

Examining the motion of a tracer particle in comparison to the motion of the centre 
of mass of the polymer chain in the bath of SPPs, we observe a giant diffusivity for the 
driven polymer. We ascribe this to the parachute like shape of the polymer in the SPP 
bath. Due to the accumulation of SPPs in the concave region of the chain the polymer 
performs an extended ballistic motion over time scales, that are considerably longer 
than that of a single monomer. The chain motion at long times becomes Brownian, 
but with an unexpectedly high diffusivity. Interestingly, even at time ranges on which 
the time averaged MSD is super diffusive, the distribution of the particle displacement 
remains Gaussian. This result differs from experimental observations of an extended 
exponential decay of the displacement of a tracer in swimming microorganisms It 
would thus be interesting to check whether incorporation of hydrodynamic interactions 
would reproduce such non-Gaussian behaviour in the model of SPP baths. Moreover, 
it is an interesting questions whether the effect of SPP pooling and the ensuing giant 
diffusivity enhancement of the polymer motion also arises in three dimensions. 

Our analysis here was performed with an in vitro bath of SPPs in mind. In 
particular the crowding fraction was chosen to be fairly low. Inside living cells active 
particles such as molecular motors drive the environment out of equilibrium and 
fluctuating forces inside cells may indeed become an order of magnitude larger than 
at the conditions of thermal equilibrium [16]. Goncurrently, the metabolic cell activity 
signihcantly affects the nature of the cytoplasm m- However, the (macromolecular) 
crowding fraction in cells typically is of the order of 30 ... 35% [IHl Il9] and thus exceeds 
the value of our simulations by far. Moreover, these crowders are quite complex 
macromolecular objects, which tune the reaction kinetics stability of biopolymers [50] . 
It will therefore be interesting to extend our study to higher crowder fractions and 
different particle geometries, such as star shapes [5T] . 
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